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The accurate determination of the thermal properties of biomass materials is particularly important 
when studying biomass combustion processes. This paper reports a method for the determination of the 
maximum sampling error and the confidence intervals of thermal properties obtained from TG-DSC 
analysis. Special care must be taken with the sampling procedure to achieve an acceptable degree of 
error and low statistical uncertainty. The obtained levels of uncertainty and error demonstrate that the 
properties evaluated by the analysis were representative of the fuels. A study of the trends and time 
correlations was performed to ensure that the observed correlations arise from completely random 
causes. The obtained results are particularly interesting for biomass energy applications. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Global concern about environmental protection grew consid¬ 
erably in the last few decades to culminate in the Kyoto Protocol [1 ], 
which set major directives and acceptable pollutant emission 
levels. The participants in the 2009 United Nations Climate Change 
Conference recognized that climate change is one of our primary 
challenges and that actions should be taken to avoid any global 
temperature increases. The Kyoto protocol, which expires in 2012, 
was ratified at the Copenhagen Summit so that its policies may be 
continued; however, no agreement was reached with respect to 
emissions reductions beyond 2012, and the final agreement is not 
legally bound. After the ratification of the Kyoto Protocol, consid¬ 
erable efforts have been made to measure and control pollutant 
output from all energy processes and especially to minimize 
greenhouse gases. The development of renewable, clean energies 
such as biomass has become an important working area as part of 
the proposed action plan [2,3], Biomass cofiring has, in fact, become 
one of the most worthwhile methods to reduce pollutant emissions 
from energy-production processes because the adjustments 
required to power plants entail low costs [4—6], 
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Several different technologies are normally applied in cofiring 
processes [7], The main advantages of each technology have been 
highlighted by various authors [4,5,7-12], Cofiring can help to 
reduce specific emissions of CO2 because of the closed carbon cycle; 
the low sulfur content of the biomass helps minimize SO2 emis¬ 
sions, and NO* emissions also show a positive trend. Cofiring 
increases the operational flexibility of the process, which reduces 
dependence on fossil fuels, such as coal. However, the main 
drawbacks of this technology are the additional cost of adapting 
combustion facilities and the increased fouling and corrosion of 
equipment [7], 

Knowledge of the chemical composition, thermal behavior and 
reactivity of biomass is essential for the effective design and 
operation of thermochemical conversion units [13]. Thermoana- 
lytical techniques, such as thermogravimetric analysis (TG) and 
differential scanning calorimetry (DSC), provide this information in 
a straightforward manner [14-18]. TG analyses are based on the 
volatilization rate of fuels, which is dependent on the heating rate 
applied to the sample and the type of fuel. DSC, in contrast, is based 
on measurements of the energy change of the sampling unit under 
a predetermined heat cycle [19,20], 

The intrinsic heterogeneity of the biomass and the small amount 
of sample used in TG/DSC experiments makes it difficult to deter¬ 
mine the thermal properties accurately. Therefore, to determine the 
characteristics of the biomass with an acceptable, clearly defined 
level of uncertainty, a well-defined TG/DSC method must be 
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developed. Many studies on the accuracy of TG/DSC experiments 
have been published [21—25], and various sampling methods have 
been proposed [26,27], Currently, TG/DSC methodologies are often 
based on small samples obtained from large batches; thus, careful 
reduction is necessary to prevent segregation and stratification 
problems [24], In general, a good sampling method should be able 
to achieve a representative sample without being affected by the 
aforementioned problems. 

A new methodology for the sampling of solid biomass, including 
the determination of error associated with the measurement of 
thermal properties, has been presented and validated [28—31]. This 
method is independent of the origin, appearance and packaging of 
the batch used to acquire the samples. 

The statistical uncertainty associated with eight different 
properties of solid fuels, which has been calculated from moisture 
and ash tests, was presented in [30], The values obtained show that 
despite the heterogeneity of the fuel itself, a well-planned 
campaign of samples can extrapolate the properties of the 
samples from the entire lot with an uncertainty controlled and 
quantified. 

A methodology for prompt analysis of heterogeneous solid fuels 
was developed taking special care with the sampling procedure to 
achieve an acceptable degree of error and low statistical uncer¬ 
tainty [31 ]. Two approaches for the propagation of errors are also 
given and some comparisons are made in order to determine which 
may be better in this context. Acceptable levels of uncertainty were 
reached demonstrating that the samples obtained in the process 
are representative of the overall fuel composition. 

The determination of the maximum sampling error and confi¬ 
dence intervals of thermal properties obtained from thermogravi- 
metric analysis (TG), including moisture, volatile matter, fixed 
carbon and ash content was studied [28,29], The results of the 
previous studies were compared to those of a prompt analysis, and 
a correlation between the mean values and maximum sampling 
errors of the methods were not observed. In general, low and 
acceptable levels of uncertainty and error were obtained with this 
methodology, demonstrating that the properties evaluated by TG 
analysis were representative of the overall fuel composition. 

In the present study, the error associated with the aforemen¬ 
tioned methodology and the confidence intervals of the specific 
heat of dry biomass and coke were determined, along with the heat 
of combustion of coke. The overall uncertainty of the measure¬ 
ments was defined, which allowed us to determine the minimum 
number of samples necessary to achieve an acceptable level of 
reliability. A comparative study on the mean values of the thermal 
properties and the corresponding uncertainties in TG experiments 
was conducted [28,29], and a relationship was not observed. 
Moreover, the confidence level and error associated with the 
measurement of thermal properties were not well correlated. 

2. Experimental section 

All materials were handled under the same conditions; they 
were tested in the same laboratory by the same analyst. Although 
the materials were exposed to environmental conditions for a short 
period of time (less than 30 min) during their manipulation, the 
effects of environmental variations in the properties of the mate¬ 
rials were not taken into account because the variations in 
temperature and relative humidity were insignificant over such 
a short time period. 

Laboratory instruments were verified and calibrated to assure 
that the experimental methodology was accurate. Errors registered 
during the experiments were considered to be random errors and 
were related to the precision of the experiment. Thus, these errors 
were quantified in the total sampling error. 


2.1. Materials 

Several different materials from agriculture and forestry were 
selected for the study. The materials cover a broad spectrum of solid 
biomass that could be used as fuel in cofiring processes. The agri¬ 
cultural materials were stored in large bags, and the forestry 
materials, in pellet form, were stored in sacks. The materials of 
agricultural origin were pine kernel shells, almond shells, hazelnut 
shells and ground olive stones. The materials of forest origin were 
pine pellets, oak pellets, brassica pellets and poplar pellets. 


2.2. Sampling and reduction of the samples 

Depending on the material, the sampled masses varied from 
320 10 -3 kg to 730-10 3 kg. Fuel samples were obtained from 
a tube sampler that was designed according to the requirements 
specified in CEN/TS [32] and the work of Pierre Gy [33—35], The 
sampling methodology used to obtain the fuel sampling units and 
the method used to reduce the samples are described elsewhere 
[28-31 ]. Fuel samples were obtained through a tube sampler that 
was designed to work with various types of solid biomass. The 
nominal maximum size, d, of the material sampled is taken as 
20 mm [31], so the tube sampler should be able to collect at least 
V min = 0.05 • d = 0.05 • 20 = 1 dm 3 = 1(T 3 m 3 [31]. Table 1 shows 
the average weight of the sampling units selected for TG/DSC 
analysis. Tweezers were used to place the sampling units into the 
crucibles. 

2.3. TG/DSC analysis methodology 

All TG/DSC experiments were performed on a TG-DTA/DSC 
SETARAM Labsys electronic thermobalance, which can achieve 
a maximum temperature of 1600 °C and heating rates from 0.001 to 
50 °C min -1 . To avoid heat- and mass-transfer limitations, 
approximately 20 • TO -6 kg of sample was used, and platinum 
crucibles without lids were employed. All experiments were 
initially conducted under a flow of nitrogen at a rate of 45 mL min -1 
to avoid oxidation of the sampling elements. After a temperature of 
600 °C was reached, dry air was introduced instead of nitrogen to 
allow an estimation of the heat of combustion of the coke that 
remained in the container. The temperature cycle and gas 
employed at each temperature are shown in Table 2. At the 
beginning of the tests, all of the equipment and sampling units 
were in thermal equilibrium at 30 °C. 

Moisture, volatile, char and ash content can be by determined 
TG/DSC results [28,29], but the main objective of this paper is to 
determine the thermal properties of the studied materials, 
including the specific heat (c p ) of the dry biomass and coke and the 
heat of combustion of the coke. 

The TG curve obtained is shown in Fig. 1, which also shows the 
DSC curve. The TG curve represents the evolution of the mass as 
a function of the temperature. The DSC curve allows the calculation 


Table 1 

Average weights of the sampling units selected for the tests. 






unit weight [kg] 


Almond shell (As) 
Ground olive stone (Gos) 
Oak pellets (Op) 

Pine Pellets (Pin) 
Hazelnut shell (Hs) 

Pine nut shell (Pns) 
Poplar pellets (Pp) 
Brassica pellets (Bp) 


21.53 x ltr 6 

22.44 x 10- 6 

21.44 x 10- 6 
20.70 x 10- 6 
21.29 x lO 6 
21.98 x 10 6 
22.46 x 10- 6 
21.18 x 10“ 6 
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Table 2 

Thermal evolution of the samples in the TG/DSC 
WTO TendTO S^Cmm" 1 ] 
30 70 30 

70 90 15 

90 105 2 

105 105 0 

105 145 10 

145 145 0 

145 500 10 

500 500 0 

500 600 20 

600 600 0 

600 600 0 

600 700 20 

700 700 0 


240 

600 

2130 

3600 

300 

600 

2400 

300 

600 


Atmosphere 

~N ~ 2 


n 2 

n 2 

n 2 

n 2 

n 2 

n 2 


n 2 

Air 


of the flow of energy by integrating the surfaces of the positive or 
negative peaks for exothermic or endothermic processes, 
respectively. 

The energy absorption is calculated from the analysis of the 
endothermic peaks that occur during the test between 105 °C and 
145 °C for the dry biomass and between 500 °C and 600 °C for the 
coke. The following Equation (1) allows the calculation of the 
specific heat (c p ) of dry biomass and coke for every presented 
material: 


d H 
dt 


d T 
C ”dt 


(1) 


The energy of coke combustion is provided by the analysis of the 
exothermic peak that occurs in the vicinity of 600 °C when dry air is 


supplied instead of nitrogen to begin the combustion process. 
Specific values are obtained by taking into account the loss of mass 
during the process. 

2.4. Statistical treatment 
2.4.1. Measured variable error 

According to Gy [33], the batch to be sampled can be considered 
as a zero-dimensional object. Using the previously described 
sampling procedure, a sample of the total batch was chosen. 
Assuming that the sampling procedure is correct, the sampling 
error, SE, can be expressed as the sum of two components: the 
fundamental error, FE, and the segregation and grouping error, SGE. 
Moreover, because these two errors are independent, the following 
relationship between their variances holds: 


o- 2 (SE) = <t 2 (FE)+<7 2 (SGE). (2) 

The FE is related to the constitutional heterogeneity; it is never 
zero, and it represents the minimum sampling error that can be 
made. The variance of the fundamental error can be expressed as: 

^ (FE) - (sL-iBtK (3 > 


where M m is the sample mass, Mi is the mass of the whole batch 
and Hli is the heterogeneity invariant given by: 




Fig. 1. The obtained DSC and TG 
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Table 3 

Values obtained from the TG/DSC analysis for c Pi db [J g _1 K -1 ], c PX k LI g _1 K -1 ] and h CO mb [kj g -1 ]. 












#1 #2 #3 

1.265 1.360 1.469 

2.345 2.373 2.018 

29.892 27.710 30.328 

1.798 1.353 1.506 

2.712 1.916 3.130 

33.494 29.690 28.991 

2.158 2.334 0.678 

2.943 2.622 

25.383 27.000 25.457 

1.680 1.715 1.612 

0.858 3.298 2.643 

32.849 32.540 31.178 

2.691 2.920 2.748 

5.565 4.494 

27.615 27.551 27.904 

1.698 1.845 1.904 

4.149 3.936 4.302 

27.717 27.971 27.439 

1.333 1.436 1.286 

2.018 1.836 1.647 

30.087 30.996 30.864 

1.532 2.435 1.381 

1.966 - 0.468 

28.950 29.033 28.850 


#4 #5 #6 

1.559 1.470 1.475 

2.468 2.370 2.138 

29.149 27.711 27.749 

1.371 1.419 

2.319 2.299 

30.810 28.236 

1.133 1.150 1.351 

0.800 

25.574 26.125 26.334 

1.630 1.801 1.697 

2.922 2.634 2.199 

31.661 30.841 31.370 

2.986 3.268 2.852 

4.299 4.990 4.067 

27.670 27.613 27.736 

1.762 1.765 1.717 

3.471 3.263 3.277 

27.153 26.768 26.497 

1.382 1.056 1.318 

1.437 3.151 2.234 

31.104 29.955 30.082 

1.368 1.412 

0.794 

27.741 28.356 


#7 #8 #9 

1.322 1.853 1.806 

2.608 2.840 

29.638 29.602 29.878 


2.232 2.241 2.378 

2.775 2.814 2.577 

26.404 26.815 27.339 

1.645 1.799 1.837 

2.185 2.158 3.093 

29.792 29.934 28.948 

2.863 2.838 2.243 

3.891 4.147 3.327 

26.707 26.990 25.570 

1.737 1.702 1.737 

2.833 2.738 3.040 

27.001 26.559 26.835 

1.387 1.413 1.413 

3.475 2.206 2.074 

31.322 31.386 30.412 


#10 


2.236 


1.155 

1.514 

26.998 


1.509 

2.395 

29.073 

1.489 

2.475 

30.244 

1.739 

2.422 

26.270 

1.713 

2.443 

31.012 

2.764 

4.348 

27.058 

1.702 

3.252 

27.094 

1.336 

2.231 

30.690 

1.076 

28.586 


S 2 


0.041 

1.121 

0.033 

0.213 

4.197 

0.428 

0.649 

0.494 

0.007 

0.523 

1.653 

0.101 

0.470 

0.901 

0.041 

0.659 

0.235 

0.013 

0.449 

0.316 

0.209 

0.621 

0.292 


where a, is the value of a certain property in a fragment or sampling 
unit, ai is the value in the whole batch, M,- is the mass of the frag¬ 
ment and N F is the number of fragments. 

In view of the above expressions, it is easy to deduce that the 
variance of the fundamental error is zero if, and only if, the sample 
is the whole batch, i.e., M m — Ml, or the material is completely 
homogeneous, i.e., a* = ai, i = 1,2,...,% 

The SGE is related to the distributional heterogeneity. The 
variance of the SGE, i.e., <j 2 (SGE), cannot be calculated. However, 
because the relationshipO < (^(SGE) < (^(FE) has been verified [36], 
then the following relationship can be easily deduced: 


M m > 7.68—^—. (9) 

SE 2 ^ 

More details about these results are available elsewhere [27], 

2.4.2. Determination of confidence intervals 

Another objective of this study was to approximate ai, which is 
the value of a property in the entire batch. We assume that as is 
similar to the value of the property in the selected sample, i.e., 
a random variable with the following distribution: 

a s ~N(a L ,a(a s )). (10) 


<7 2 (FE) < <t 2 (SE) < 2<t 2 (FE). 


(5) 


Assuming that the sampling error follows a normal distribution, 
i.e., SE ~ N(0 ,<j(SE)), we can ensure with a confidence level of 95% 
that: 


(SE) < 1.96n/2(t(FE) = 1.96^2^ (^- - T-)h/ l . 


(6) 


Finally, assuming that M m < Mr,: 


From the definition of the SE and Equation (5), an approxima¬ 
tion of the variance was obtained: 

<7 2 (a s ) = <r 2 (SE)-a 2 < (11) 

Finally, the value of a property in the entire batch, ai, can be 
estimated as the mean of the experimental values, i.e., a realization 
of a s , and confidence intervals for a L can be obtained from the 
approximation of the variance of as. 


|SE| < 1.96 


j2 ■ HIi 

V M^T 


(7) 


Using the latest inequality, some useful bounds, with a confi¬ 
dence level of 95%, can be inferred for the SE and the mass of the 
sample: 


3. Results and discussion 

Values obtained during the TG/DSC tests of the different mate¬ 
rials to obtain the specific heat of wood (c p ,db). the specific heat of 

Table 4 


1 If the mass of the sample is constant, the SE exhibits an upper 
bound of the maximum sampling error given by: 


|SE| < SEmax 



(8) 


2 If a maximum sampling error is set, then the mass of the 
sample should be: 


Intrinsic heterogeneity of c p> db, c PX k and h CO mb of different biomass materials. 


*-p,db 

~As 1.61 x 10 2 

Bp 1.20 x 10“ 2 

Gos 1.26 x 10 _1 

Hs 2.03 x 10- 3 

Op 1.19 x 10“ 2 

Pin 1.28 x 10“ 2 

Pns 6.72 x 10“ 3 

Pp 6.32 x 10~ 2 


1.00 x 10" 2 
2.79 x 10~ 2 
9.48 x 10“ 2 
7.78 x 10- 2 
2.24 x 10- 2 
5.61 x 10- 2 
8.02 x 10“ 2 
3.57 x 10 _1 


1.18 x 10“ 3 
3.67 x 10" 3 
6.36 x 10" 4 
1.53 x 10“ 3 

2.88 x 10- 4 
2.99 x 10“ 4 
2.86 x 10~ 4 
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Table 5 

The minimum sample mass (expressed as number of sampling units) required to achieve a predetermined maximum sampling error for the determination of cp.db. 


Minimum sample size required to obtain a determined maximum error 


Maximum error 


HI L 

0.001 

0.005 

0.01 

0.05 


As Bp Gos Hs Op Pin Pns 


1.61 x 10“ 2 
1.24 x 10 5 
4.95 x 10 3 
1.24 x 10 3 
4.95 x 10 1 


1.20 x 10“ 2 
9.23 x 10 4 
3.69 x 10 3 
9.23 x 10 2 


1.26 x 10-’ 
9.66 x 10 5 
3.86 x 10 4 
9.66 x 10 3 
3.86 x 10 2 


2.03 x 10“ 3 
1.56 x 10 4 
6.24 x 10 2 
1.56 x 10 2 
6.24 x 10° 


1.19 x 10“ 2 
9.12 x 10 4 
3.65 x 10 3 
9.12 x 10 2 
3.65 x 10 1 


1.28 x 10“ 2 
9.82 x 10 4 
3.93 x 10 3 
9.82 x 10 2 
3.93 x 10 1 


6.72 x 10“ 3 
5.16 x 10 4 
2.06 x 10 3 
5.16 x 10 2 
2.06 x 10 1 


_PP_ 

6.32 x 10“ 2 
4.85 x 10 5 
1.94 x 10 4 
4.85 x 10 3 
1.94 x 10 2 


Table 6 

The maximum sampling error SE nax that corresponds to a given sample mass (expressed as number 


Sample size 


HIi 


200 


Maximum error 



as a function of a given sample size 


_Bp_ 

1.20 x 10“ 2 
3.04 x 10-' 
9.61 x 10“ 2 
3.04 x 10“ 2 
2.15 x 10“ 2 


Cos 

1.26 x 10-’ 
9.83 x 10“’ 
3.11 x 10 1 
9.83 x 10“ 2 
6.95 x 10“ 2 


Hs 

2.03 x 10“ 3 
1.25 x 10“’ 
3.95 x 10 2 
1.25 x 10“ 2 
8.83 x 10“ 3 


of sampling units) for the determination of c p ^b. 


Op Pin 

1.19 X ltr 2 1.28 X 10“ 2 

3.02 x 10“’ 3.13 x 10“’ 

9.55 x 10“ 2 9.91 x 10 2 

3.02 x 10“ 2 3.13 x 10“ 2 

2.14 x 10“ 2 2.22 x 10“ 2 


Pns 

6.72 x 10“ 3 
2.27 x 10“’ 
7.18 x 10“ 2 
2.27 x 10“ 2 
1.61 x 10“ 2 


_PP_ 

6.32 x 10“ 2 
6.97 x 10- 1 
2.20 x 10 1 
6.97 x 10 2 
4.93 x IQ- 2 


coke (Cp.ck) and the heat of reaction (h C omb) are shown in Table 3, 
which also includes the mean and variance of the variables. 

Table 4 presents the values of HIi, that were calculated according 
to the method described in section 2. For calculations included in 
this table, fragment is assumed to be a dimensionless unit of mass, 
Mi = 1, so that the mass sample is represented as NF sampling units. 
That is to say, for each material the unit of mass stands for the 
sampling unit shown in Table 1. 

Tables for any of the properties studied can be generated from 
the heterogeneity invariant so that the maximum error of a sample 
can be determined for a known sample size. The minimum sample 
size necessary to obtain the previously set maximum error of 
a sample can also be calculated. The minimum sample sizes 
necessary for a given maximum error and the maximum errors for 
a given sample size are presented in Table 5 and Table 6, respec¬ 
tively. The utility of these tables and examples of how they can be 
used are discussed elsewhere [28—31 ]. 

Following the methodology described in section 2.4, confidence 
intervals with a 95% confidence level for c Pi db, c p , c k and /i C omb of each 
material were calculated. Table 7 shows the limits of these confi¬ 
dence intervals. The ranges are generally quite small, i.e., on the 
order of tenths, except for the /t C omb of brassica pellets (Bp), the 
hcomb of hazelnut shells (Hs) and the c PfC k of poplar pellets (Pp), 
which exhibit ranges on the order of units. 

An analysis of Table 3 reveals that it was not possible to deter¬ 
mine the studied properties for all of the tested samples. The c PiC k of 
Sample #7 for almond shells (As), Samples #3, #5 and #6 for 
ground olive stones (Gos), Samples #1 and #10 for oak pellets (Op) 
and Samples #2 and #5 for poplar pellets (Pp) could not be 
determined. A notable difference was also observed between the 
minimum and maximum variance of the properties studied. 


To relate the dissimilar values of the variance associated with 
average values that vary widely from each other, it is appropriate to 
compare the values of the coefficients of variation that allow an 
assessment of the magnitude of the dispersion from the mean. The 
coefficients of variation for each property of each of the materials 
are listed in Table 8. 

Coefficients of variation (Table 8) for the specific heat of dry 
biomass and for the specific heat of coke are on the order of tenths. 
The specific heats of dry biomass are significantly lower, except for 
in the cases of almond shells (As) and ground olive stones (Gos). 
The coefficient of variation of the specific heat of coke for pine nut 
shells (Pns) reaches a value below the hundredths. 

Because the measurements of the properties may exhibit 
temporary features, an additional analysis was performed to check 
for an underlying time-series structure. Because the TG/DSC tests 
were performed for each sampling unit, the derivation from each 
test of the values of the different properties, i.e., moisture, volatile, 
fixed carbon and ash, which were used in a previous work, and the 
specific heat of dry biomass, specific heat of coke and heat of reaction 
of coke, allows a comparison of the linear trend and random varia¬ 
tion of the seven properties with previously published results [28], 

Table 9 summarizes the results of statistical analyses applied to 
the sample data, i.e., Pearson correlation tests for the presence of 
a linear trend and the Ljung—Box test for the correlation time of 
random variation. The first two columns are used to verify the 
existence of a linear trend using the Pearson correlation coefficient 
and provide the corresponding p-value, respectively. The third and 
fourth columns show the p-values of the Ljung—Box test, which is 
a statistical hypothesis test used to test the null hypothesis that the 
residuals of a time series are uncorrelated. 


Confidence intervals < 


properties. 


_ Wl g- 1 K- 1 ] 

As 1.51 ± 1.7710-' 

Bp 1.49 ±2.02 10“’ 

Gos 1.74 ± 5.70*10“’ 

Hs 1.71 ± 7.13-10- 2 

Op 2.76 ±2.64 10“’ 

Pns 1.34 ± 1.0110“’ 

Pp 1.63 ± 5.07 *10“’ 


Cp.ckLI g“’ K- 1 ] 

2.39 ±2.35*10“’ 
2.48 ± 5.12*10“’ 
2.42 ±7.81*10“’ 
2.44 ±6.30*10“’ 
435 ±5.7*10“’ 
3.25±6.75*i0“’ 
2.23 ±5.84*10“’ 


frcomblH g^] 

29.07 ±9.22-10-’ 
30.24 ± 2.27 
26.27 ±6.1210“’ 
31.01 ±1.12 
27.06 ±7.89 10“’ 
27.09 ±4.03 10“’ 
30.69 ±4.90 10“’ 
28.59 ±5.9910“’ 


Table 8 

Coefficients 


of variation. 


Properties 


As 0.135 0.107 

Bp 0.122 0.186 

Gos 0.376 0.333 

Hs 0.049 0.296 

Op 0.115 0.158 

Pin 0.119 0.250 

Pns 0.085 0.300 

Pp 0.281 0.732 


0.036 

0.068 

0.088 

0.041 

0.035 

0.018 

0.018 

0.019 
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0.147 



0.317 

0.209 

0.349 


Pp 


For a significance level of a = 0.05, the values for oak pellets (Op) 
and pine pellets (Pin) exhibit a significant trend for the specific heat 
of coke and the heat of reaction. The values of the dry biomass 
specific heat for the almond shells (As) and the heat of reaction for 
hazelnut shells (Hs) also appears to exhibit a significant trend. 
These correlations seem to be random, because they appear in 6 of 
24 cases studied, all with a limited number of tests. 

After the trend has been removed, the Ljung-Box test detects 
the correlation in time. In the case of the heat of reaction for the 
pine nut shells, when considering two lags, we conclude that this 
relationship is irrelevant because the two lags only appear in 
a property of a material. The significance of this correlation in 
a single case suggests that this result is uncertain. 

4. Conclusions 

In this paper, samples obtained following the methodology 
described in [28] were investigated by TG/DSC analysis. The heat of 
reaction of the coke, the specific heat of the dry biomass and the 
specific heat of the coke were determined for some solid biofuels. 
Statistical analyses and other tests used in previous reports [29] to 
study the linear trend and the random variation of the properties 
were employed to analyze the results obtained during the 
experiments. 

The lack of values for the specific heat of coke in samples of 
various materials indicates that the analysis has certain weaknesses 
and that the data are of limited reliability. This conclusion is sup¬ 
ported by the variances and some intrinsic heterogeneities that 
lead to wide confidence intervals and high coefficients of variation 
in virtually all cases. As a result of the confidence intervals of the 
specific heats and the heats of reaction, we estimate that the 
average values can be employed in all cases, except for the specific 
heat of coke for pine pellets. 

When comparing the linear trends and random variations of the 
properties of the eight materials listed in Table 9 with the results of 
the same test for other properties published elsewhere (Table 14 in 
previous work [28]), we observe that the materials with properties 


this tendency in all cases. Therefore, these comparisons confirm that 
the correlations that appear are derived from random causes. The 
TG/DSC analyses presented in this paper and in previous works, as 
well as the process of extraction, reduction and sample testing, show 
neither linear trends nor correlations of time. This latter conclusion 
confirms the initial supposition of zero-dimensional samples. 
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